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O ; ABSTRACT 

. We present high-resolution hydrodynamical simulations of the degree and direction of 

polarisation imprinted on the CMB by the Sunyaev-Zel'dovich effect in the the line of 
I sight to massive galaxy clusters. We focus on two contributions to the electron rest- 

frame radiation quadrupole anisotropy in addition to the intrinsic CMB quadrupole 
that contributes most of the induced CMB polarisation: the radiation quadrupole seen 
by electrons due to their own velocity in the plane normal to the line of sight, and the 
^y-^ ' radiation quadrupole due to the thermal Sunyaev-Zel'dovich effect, which is generated 

Q>^ I by a previous scattering elsewhere in the cores of the local and nearby clusters. We 

psj , show that inside the virial radius of a massive cluster, this latter effect, although being 

' second order in the optical depth, can reach the nK level of the former effect. These two 

I effects can, respectively, constrain the projected tangential velocity and inner density 

CO . profile of the gas, if they can be separated with multi-frequency observations. As the 

' information on the direction and magnitude of the tangential velocity of the ICM 

r~| . gas is combined with the temperature Sunyaev-Zel'dovich effect probing the projected 

^ i' radial velocity it will be possible to improve our understanding of the dynamics of 

the ICM. In particular, future polarisation observations may be able to trace out the 
filamentary structure of the cooler and less dense ionized gas that has not yet been 
^ incorporated into massive clusters. We also discuss how the collapse of a cluster should 

^ I produce a peculiar ring-like polarisation pattern which would open a new window to 

detect proto-clusters at high redshift. 
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1 INTRODUCTION 

Recent measurements of the TE power spectrum of the in- 
trinsic polarisation of the CMB by WMAP have put strong 
constraints on the total optical depth r to the last scattering 
surface and on the possible nature of the sources responsi- 
ble for reionisation. These measurements have also provided 
impetus for the further development of polarisation-sensitive 
receivers and for the theoretical modelling of the foregrounds 
contaminating the CMB polarisation. 

On the other hand, new generation X-ray satellites such 
as Chandra or XMM-Newton have produced a wealth of high 
resolution maps of the intracluster medium (ICM) of nearby 
clusters. They support some aspects of our previous picture 
of the hydrostatic state of the plasma, but they invalidate 
others, and, most interestingly, bring out new phenomena 
such as cold fronts. However, it is well known that detailed 
studies of the cluster ICM in the X-ray band constrain the 
innermost parts, while the dynamics of the gas is also in- 
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teresting at larger radii where for example gas flowing in 
from filaments, pristine gas accreting from the low-density 
regions of the Universe, or subclumps of cold gas of merging 
units encounter the hot ICM. 

The Sunyaev-Zel'dovich effect (hereafter, SZ) encapsu- 
lates the processes of Thomson and inverse Compton scat- 
tering of CMB photons by the free electrons of the ICM (see, 
e.g. Sunyaev & Zeldovich 1980a; Carlstrom et al. 2002). In 
CMB temperature maps, both the thermal SZ (hereafter 
tSZ) effect and the kinetic SZ (hereafter kSZ) effects con- 
tribute to the secondary anisotropies. With CMB tempera- 
ture variations depending at fixed observing frequency only 
on the gas pressure integrated along the line-of-sight, the 
tSZ effect has long been recognized as a powerful probe for 
studying the morphology of the ICM out to a substantial 
fraction of the virial radius of the cluster. The kSZ effect 
typically amounts to a few percent of the thermal effect and 
maps the bulk velocity of the ICM projected along the line 
of sight, to the extent that the ICM can be considered as 
a coherent bulk. Multi-frequency observations are necessary 
to separate the kSZ effect from the tSZ effect and extract 
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information on the radial velocity of the gas. Recently, us- 
ing Eulerian hydrodynamical simulations, Nagai & Kravtsov 
(2003) have studied the bias due to the internal motion of 
the gas in a massive, 2.4 x lO^^'/i"^ Mq cluster on the es- 
timates of the cluster bulk velocity along the same lines of 
sight where one measures the kSZ effect. While they focused 
on the ability to determine the global motions of massive 
haloes to later constrain cosmological parameters by means 
of the cluster pairwise peculiar velocities, it is also possible 
to constrain the motion of the ICM in the rest frame of the 
cluster. In fact, the kSZ effect and some of the SZ polarisa- 
tion signals are complementary as they allow one to obtain 
information on, respectively, the projection along the line of 
sight of the density-weighted radial velocity of the gas, and 
the projection (amplitude and direction) along the line of 
sight of the density-weighted tangential velocity of the gas. 

On CMB polarisation maps, the SZ-induced polarisa- 
tion comes mainly from four terms. The first, which we want 
to emphasize, comes from single scattering of the quadrupole 
generated as radiation is Doppler-shifted due to the tan- 
gential motion of the gas (hereafter kpSZ), the second is 
due to single scattering of the CMB intrinsic quadrupole 
(hereafter ipCMB), and the last two are second scatterings 
of quadrupole anisotropics induced in the CMB by respec- 
tively previous t- and kSZ effects induced by first scatterings 
(hereafter respectively p^tSZ and p'^kSZ ), when the y pa- 
rameter/optical depth seen by the last scattering electron is 
not isotropic. These last two effects arc potentially powerful 
probes of the central g£is profiles in clusters. 

For last scattering electrons inside and outside clusters, 
the p^kSZ effect is negligible relative to the p'^tSZ effect if 
clusters are sufficiently hot and massive: these terms scale 
respectively as ?7T^ and /3r^ (Sazonov & Sunyacv 1999, here- 
after S99) where = kT/mc^ and (3 the tangential veloc- 
ity. Indeed, massive clusters have typically 77 ~ 3 10~^ and 
/3 ~ 3 10"'\ While the ipCMB can reach levels of order 50 
riK (S99), it is possible to eliminate it either by restricting 
oneself to local clusters at z < 0.3 located in the directions 
of the z=0 CMB quadrupole where the ipCMB vanishes, or 
by using intensity maps of the tSZ efi^ect and a model for 
the temperature distribution of the ICM. We will not tackle 
this term in detail as we consider that it is well understood 
and straightforward to simulate. 

However, the relative amplitude of the ktSZ and the 
p^tSZ effects inside and in the vicinity of massive clusters 
calls for a more quantitative assessment as this is expected 
to depend strongly on the density and on the temperature 
profiles of the cluster. We will analyse the p^tSZ efl^ect at 
some length, for three reasons. Firstly, the degree of CMB 
polarisation due to this second scattering of the tSZ has 
only been evaluated analytically, using spherically symmet- 
ric models for clusters (S99). Secondly, as pointed out by 
Sunyaev & Zeldovich (1980b, hereafter S80), it is interest- 
ing as the second scattering of the tSZ effect depends only on 
the local free electron density and not on the temperature, 
similarly to the scattering of the intrinsic CMB quadrupole. 
It can therefore be viewed as a means of detecting large- 
scale structures other than clusters: groups and filaments 
where the inner gas density is a few times the mean cosmic 
value. In this regime, its key advantage over the polarisation 
induced by the intrinsic CMB anisotropy resides in that it 
can be a factor of 10 more eSicient, depending on observing 



frequency. Thirdly, and most importantly, double scattering 
can provide information about the inner gas density profiles 
of clusters, cross-checking the constraints derived from di- 
rect X-ray observations, which depend on a model for the 
local temperature of the gas when it is not possible to mea- 
sure it directly. 

Our hydrodynamical simulations show that maps of 
CMB polaxisation can provide useful information about the 
profiles and dynamics of the ICM gas and constrain our 
picture of gas infall, mixture and heating as well as the dy- 
namics of cold subclumps. They indicate a possible way to 
detect proto-clusters in their early stage of collapse. Besides, 
they show how the cluster SZ thermal effect can polarise the 
CMB over the filaments of gas between clusters. 

The layout of this paper is as follows. In section 2, we 
review the analytical formulae we employ to compute the 
degree and direction of polarisation, due to kpSZ effect and 
to second scatterings of the tSZ effect, i.e. the p'^tSZ effect. 
We present a semi-analytical description of the latter which 
enables us to obtain a substantial gain in computational 
time. We also summarise our assumptions. In section 3, we 
describe the setup of the simulations and the construction of 
the polarisation maps. Section 4 presents the strongest de- 
pendences of our results on the cluster gas profile. Section 5 
discusses a general observational methodology to deal with 
contamination efi^ccts in order to retrieve the original pro- 
jected density-weighted tangential and radial velocity field. 
Conclusions are presented in section 6. 



2 ANALYTICAL DESCRIPTION 

In this section, we briefly summarise the analytical expres- 
sions we use for the main terms contributing to the polar- 
isation of the CMB induced by the SZ effect discussed by 
S80 (see also S99, Chluba & Mannheim 2002). 

We fix an ortlionormal frame (X, Y, Z ) at the point P 
where the polarisation is induced: it is the position of the 
scattering electron in the case of single scattering, or the 
position of the last scattering electron in the case of double 
scattering. We choose the Z direction along the line of sight, 
from P to the observer. X and Y label the sky coordinates 
(we assume an euclidian patch) . We use the Stokes parame- 
ters (/, Q, U , V) to describe the amount of polarisation as 
seen by an observer looking in a given direction. / gives the 
total intensity of the light, while Q and U measure the degree 
of linear polarisation, and V the degree of circular polari- 
sation. Thomson scattering of anisotropic but unpolarised 
incoming radiation cannot generate circularly polarised ra- 
diation so that V = Q. The degree of polarisation is given 



by: 



P- 



(1) 



and the direction of linear polarisation (/), measured with 
respect to the X axis, is: 

IJ_ 

Q' 

In the rest of this section we deal with local effects on 
the line of sight. This approach is valid since the Stokes pa- 
rameters are additive along the line of sight and the resulting 



tan{2(f>) 



(2) 
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Q and U seen by the observer read as follows: 

Q = I dQ 
U = J dU, 

where dQ and dU are local Stokes parameters at position 
(X, Y, Z ). As we deal with small optical depths throughout 
the paper we assume that the background intensity 7o of the 
CMB is not affected by SZ scattering: we set 7tot = I — lo 
to compute the degree of polarisation. 

We now present in turn (1) the polarisation induced by 
the transverse motion of the gas with respect to the line of 
sight, (2) the effects of double scattering, and (3) the con- 
tribution due by scattering the CMB intrinsic quadrupolo 
which together generate most of that part of the CMB po- 
larisation due to the SZ effect. We will mention in the last 
paragraph of section 5 two other terms which can contribute 
a large fraction of the total polarisation of the CMB at a 
given frequency: the intrinsic polarisation imprinted at last 
scattering and the main types of foregrounds, i.e. dust grains 
and point sources emitting synchrotron radiation. 



2.1 Kinetic polarisation 

The kinetic polarisation kpSZ comes from single scattering 
of a CMB photon by a free electron moving in the plasma. 
Here, the electron rest frame anisotropy of the radiation 
necessary for Thomson scattering to generate polarisation 
is entirely due to the electron velocity (S80). We refer the 
reader to S99 and S80 for the detailed derivation and we 
just summarise their results here. 

The Q and U components are given in the rest frame 
of the observer (all quantities depend on (X, Y, Z), though 
not explicitly written): 



dQ = -0.1/o.o-T/(a;.)ne/3|cos(2x) dZ 
dU = -Q.lIo^aTf{x^)nePtsm{2x) dZ 



(4) 



Here, ctt is the Thomson cross-section, rie is the local density 
of electrons, is the dimensionless parameter of the Planck 
law: = and: 



e--(e-- + l) , 



(5) 



is the spectral shape of the polarisation induced by the CMB 
in a non-relativistic approximation and simply is a transla- 
tion of the black-body spectrum because of Thomson scat- 
tering, /i — J-,^ , where 1? is the velocity of the electron and 
Pt = — M^) is its tangential component, and x is the 
angle between the tangential velocity and the X-axis in the 
tangential plane. This effect is simple to implement as one 
simply needs to know the local state of the gas: given the 
value of each field everywhere in the simulation the summa- 
tion is straightforward. 



2.2 Double scattering due to finite optical depth 

S80 noted that since the optical depth t along the line-of- 
sight through massive clusters can reach ~ 0.03, double scat- 
tering cannot be neglected. In this case, CMB polarisation 
is due to second scattering of CMB secondary anisotropics 
generated by a first scattering in a dense region, either in the 



same structure or in a nearby one. The whole process can be 
alternatively viewed as scattering of the tSZ and kSZ effects. 
However, due to the kSZ effect being typically more than an 
order of magnitude smaller than the tSZ effect for massive 
clusters ( ^ 6 x 10^'* M©), and as we will show that scatter- 
ing of the tSZ (p^tSZ ) gives a degree of polarisation of the 
same order as the kpSZ effect of the previous paragraph, it 
is safe to neglect scattering of the kSZ (p^kSZ ). S99 give an 
explicit formal formulation of this second scattering effect: 



dQ = 

dU = 

where: 



dZ^aTne(X,Y, Z) 

x / dO, sm'^{e)cos{2<j))Ah(X,Y,Z,d,(l>) 

dZ^aTne(X,Y, Z) 

xJdQ sin2(6l)sin(2^6)A7,(X, Y,'Z,6',0), 



(6) 



A7.(X, Y, Z) =yeff(X,Y, Z ,e,<t>)fT{xu)Iou, (7) 

dn = sm{e) de dcf), (s) 

frixu) is the frequency dependence of the tSZ effect: 



frixi,) = Xu 



+ 1 



(9) 



and lov is the intensity of the CMB at frequency v. and 
(j> give the direction of the radiation incoming to the second 
scattering electron located at (X, Y, Z). With our axes, 6 
is the angle with respect to Z (colatitude) while 4> is the 
counterclockwise angle to X in the (X,Y) plane, j/efi is the 
integral of rydr as seen from (X, Y, Z) in the outward di- 
rection (6, (j>) and is just the Compton y parameter in that 
direction. Interestingly, S80 pointed out that the degree of 
polarisation in double scattering at fixed tSZ effect "back- 
ground" only depends on the local optical depth, and could 
be used to probe cooler regions than possible with temper- 
ature maps of the tSZ effect itself, provided of course that 
there is already sufficient sensitivity to measure the p^tSZ 
effect in the core regions of galaxy clusters. 

In the following, we describe our computation in some 
length as we show that the p'^tSZ effect depends sensitively 
on the detailed gas density profile of clusters, that it can 
in fact dominate in amplitude over the kpSZ effect and so 
should not be neglected in single frequency observations, 
and that it could also be of interest as a means of inducing 
large-scale polarised CMB radiation. 

2.2.1 A spherical isothermal model for first scatterings 

The double integrals in (6) are difficult to compute numer- 
ically: one needs to first evaluate from each point of the 
simulation volume am integral in all directions, and then to 
perform the summation along the line of sight. Therefore, 
we have adopted a semi-analytical model: we first choose a 
spherical, isothermal model for the clusters - here a /J-model 
- to be able to compute j/e//(X, Y, Z,0, 4>) for each cluster. 
We define the /3-modeI as: 



p(r) = p* 



1 + 



-3/3/2 



(10) 
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(see, e.g., Komatsu & Seljak 2001) where Tc is the core radius 
of the cluster and p* is a normalisation constant and we take 
13 = 2/3. 

The important task is to estimate p* and Vc for each 
cluster of the simulation. The optical depth as seen from 
(X, Y, Z) in the direction {9, (j)) can be computed as: 



r(X,Y, Z,< 



= (Jt 



1 , 

P 

mpH Jo 



1 + 



dl 



1 * 2 

= ar p 

rupp, 



+ tan 



ff(X,Y, Z,g,0) 
/r,H^?-5'(X,Y, Z,e, 




(11) 



where Z is the distance from the second scattering electron 
to the point taken along the line of sight {0,4>), ro is the 
distance between the point (X, Y, Z) and the cluster lo- 
cated at (Xc, Yc, Zc), rup p is the mean particle mass of 
the ICM, and: 

fl(X, Y, Z , e, = (X - Xe) sm(9) cos(0) + 

(Y - Yc) sin(6l) sin((^) + (Z ~ Zc) cos{e). (12) 

j/eff(X, Y, Z,6,4>) is here the product of t(X, Y, Z,8,(f>) 
with the mean relative temperature r? = , where Te is 
the mean temperature of the electrons, rUe is the electron 
mass and fee is Boltzmann's constant. Once we have cho- 
sen the sot of clusters which make the most contribution 
to the tSZ effect over most of the simulation, at each point 
(X, Y, Z) we can add up the local Stokes parameters ob- 
tained from (6) and (11) for each of these clusters to obtain 
the local dQ, dU due to the p'^tSZ effect. 

The underlying assumptions here are first, the validity 
of the /3-model for the radial profile of the gas, i.e. the va- 
lidity of spherical symmetry and isothermality for the gas 
given that the X-ray surface brightness would bo fitted by a 
/9-profile, second, the validity of isothermality in the compu- 
tation of the tSZ effect itself and third, the fact that lower 
mass groups and filaments will contribute very little to the 
tSZ effect as observed from any point of the simulation. 



2.2.2 Estimating the core radius 

Because we find the relative amplitude of the p'^tSZ to the 
kpSZ effect to sensitively depend on the gas distribution, for 
massive clusters intersecting along the line of sight, we have 
computed the normalisation of the /3-model in two ways. 
This ensures that we are not affected by resolution issues 
inherent to SPH such as gravitational softening. We first 
estimate the parameters numerically under the assumption 
that the cluster gas exactly follows a /3-model, and we then 
derive an empirical best-fit to the actual simulated density 
profile of each cluster. This comparison also checks the va- 
lidity of the /3-profile for the simulated clusters. 

Numerical approach To estimate the free parameters p* 
and rc numerically, we first identify clusters using a friends- 
of-friends groupfinder (Davis et al. 1985) with linking length 



0.2 times the mean intorparticle separation. We then com- 
pute the virial mass and radius rrtvir and rvir starting from 
the most bound gas particle of the cluster. Note that we 
use only gas particles to estimate rrivir and Tvir. We define 
the virial radius as the point whore the enclosed gas density 
falls to Ac Pc /bar where Ac is the virialization overdensity 
(Bryan & Norman 1998), pc is the critical density of the 
universe and /bar = 0.13 is the mean cosmological gas frac- 
tion. We suppose that the cluster only extends up to the 
virial radius. Starting with an analytical expression for the 
virial mass myir as a function of p* and rc, it is possible to 
express the mass inside the core radius mdrc) as a func- 
tion of TTivir and rc. At fixed mvir, we then iterate to find 
the core radius: it is obtained by choosing a test radius rt, 
computing the mass m(< rt) inside rt and comparing it to 
mc{rt). Practically, we start from the most bound particle 
of the cluster and increase rt, up to the desired core radius 
where m(rt) = mc{rt) where we then set rc = rt. 



Best-fit approach As a second approach, we directly fit 
the /3-model to the gas density profile measured in the simu- 
lation. We employ an adaptive number of points depending 
on the number of particles in the cluster so that noise fluc- 
tuations on small haloes can be reduced. We always take 
more than 300 particles per halo to give a good fit, with the 
number of fitting points less than 100, and we limit ourselves 
again to the virial radius. We find reasonable agreement be- 
tween the two approaches, showing that the /3-profile and 
spherical symmetry are fair assumptions. When computing 
the p^tSZ effect in the simulations, we employ the second 
approach which gives fast results and depends on both the 
inner and outer parts of the cluster. For the massive cluster 
that we find in our high resolution simulation we compute 
mvir = 1.5 X IO^^Mq and the numerical and best-fit ap- 
proaches give respectively rc = 250 kpc and 370 kpc. 



2.2.3 Semi- analytical results 

Having set the gas density profile for the integration of (11) 

over the path of the incident radiation on the second scat- 
tering photon, we proceed to integrate (6) over (9, (f>) nu- 
merically. Figure 1 shows the results of the computation of 
Q and U, for a spherically symmetric gas distribution fol- 
lowing a /3-profile with core radius 150 kpc. Here and in all 
subsequent figures we take .x^ = 3: wo discuss the frequency 
dependence of the effects in section 4. 

As expected, the polarisation has circular symmetry. 
The signal vanishes at the centre of the cluster because of the 
particular symmetry involved at that point: the centre of the 
cluster docs not see any changes duo to any rotations of the 
space around it, so that there can be no polarisation induced 
at this point. As a consequence, the intensity pattern shows 
an annulus: we checked that the radial distance rmax to the 
circle on the annulus where the degree of p'^tSZ is maximum 
scales as rmax ~ 1.3 rc. Note that this feature has already 
been found by S99 (see their figure 3) who also integrate 
their equations numerically: their King profile with n=l is 
similar to the /3-profile assumed here. 
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Figure 1. The top images gives the value of Q and U (left and 
right panels respectively) of the polarisation of the CMB due to 
second scattering of the in-situ tSZ effect for a 150 kpc core radius 
cluster put at the centre. The figure illustrates the spatial vari- 
ations of U and Q but the real level of polarisation will include 
an overall factor dependent on the normalisation of the matter 
density profile. The gray scale is linear in Q, U and the polarisa- 
tion. Black is for the minimum relative value and white for the 
maximum relative value. The change in the intensity of the po- 
larisation is clear at the centre. As expected the axis of Q and U 
are rotated by 7r/4 and the axis of Q is orthogonal to the a;-axis. 
The global degree of the polarisation, as given hy \/ Q"^ + for 
/ = 1, is shown in the lower image. Note the annulus of high in- 
tensity around the core and the fall-off at the centre as spherical 
symmetry washes out the signal. 



2.3 Polarisation due to the CMB quadrupole 

Thomson scattering of the intrinsic quadrupole of the CMB 
by the free electrons of the ICM will contribute to the de- 
gree of polarisation observed toward massive clusters. S99 
show that this term can reach 2 x 10"'^ r ~ 4 x 10~* for 
r = 1/50, comparable and even higher than the polarisation 
due to the kpSZ or p^tSZ effects presented above. However, 
at low redshift, there are four orthogonal directions in the 
sky where one expects the degree of CMB polarisation due 
to the scattering of the quadrupole component to vanish 
(see Fig.l of S99). In practice, the directions towards which 



the CMB quadrupole vanishes have been determined with 
some accuracy by WMAP: one is near {l,b) ~ (—80°, 60°), 
a direction close to the Virgo cluster (Tegmark et al. 2003). 
Note that this only holds for local clusters {z < 0.3), as ob- 
jects at higher redshift will observe a different realisation of 
the primordial quadrupole (Kamionkowski & Loeb 1997). 



3 RENDERING A GALAXY CLUSTER 
3.1 Parameters and setup 

We use the public version of GADGET (Springel et al. 2001) to 
run the simulations after modifications to include the "stan- 
dard entropy" formulation (Springel & Herncjuist 2002), 
cooling of the gas and photoionisation by a uniform UV 
background. The simulations follow a cube of 70 Mpc co- 
moving side and assume a concordance ACDM cosmology 
with 0,0 = 0.3, h'^ = 0.019, A = 0.7, h = 0.7. We have 
not switched on cooling and photoionisation in this study 
but we will analyse their impact in a forthcoming paper ; 
here the evolution is fully adiabatic which suffices for our 
purpose. 

We constrain the density field using the implementation 
of the HofTman-Ribak algorithm proposed by van de Wey- 
gaert & Bertschinger (1996) to obtain a massive cluster at 
the end the simulation. To obtain a cluster of approximately 
1.5 X lO^'^A/0 we constrain the initial conditions to have a 
gaussian peak with a = 13.5 Mpc in the centre of the box. 
The height of the peak is three times the rms value obtained 
by smoothing an unconstrained density field with the same 
gaussian kernel. At z—Q, the optical depth of the central 
region of the cluster is ^ 3 10~^ and the mean electron tem- 
perature in the virialised gas is ~ 15 keV: the cluster is more 
massive than Coma, and probably closer in mass to Perseus. 

To speed up the simulation we have used mesh refine- 
ment when constructing the initial conditions: after run- 
ning a low resolution simulation the particles of the cluster 
have been traced back to the initial redshift to find the Lar- 
grangian region containing all the mass virialised by z=0. We 
have used two simulations: in the first the initial density field 
is realised on a 64"^ mesh embedded in a 32^ mesh covering 
the whole volume, while the second uses a 128"' mesh embed- 
ded in a 64^ mesh. In the two simulations, we have used the 
same phases for the common density perturbation modes so 
as to verify the stability of the method with respect to mesh 
resolution. We put high mass DM and gas particles on the 
nodes of the low-resolution grid and low mass DM and gas 
particles on the nodes of the high resolution grid. Following 
common practice, we set the softening length of both low 
and high resolution particles (we use the subscripts "Ir" and 
"hr" in the following) to one tenth of the their mean inter- 
particle separation. The softening length IsoU is kept fixed in 
comoving coordinates: Zsoftjr = 400kpc and Zsoft,hr ~ 30kpc. 
We have checked that the overall intensity of the efi'ects and 
the gas distribution is the same within a factor of two be- 
tween the simulations, and in the rest of the paper we will 
only discuss the simulation with high mesh resolution. 

The z=0 map of the density projected along the Z di- 
rection of the whole box is shown in Figure 2. The box is 
70 Mpc wide. The friends-of-friends groupfinder finds a lit- 
tle more than 3.5 x 10^ DM particles in the central clus- 



6 G. Lavaux, J. M. Diego, H. Mathis, J. Silk 

2 22 2.4 2.6 18 3 3.2 3.4 3.6 3.B 4 4J 4.4 4.6 4.6 5 




10 20 30 40 50 60 70 



(x) Mpc 

Figure 2. Map of the gas density of the whole high resolution 
simulation projected along the line of sight. The side is 70 Mpc 
and the cluster is located close to the middle of the box (as it 
is along the line of sight). There arc four filaments rea<;hing the 
cluster and extending at close to right angles. The colour scale is 
in decimal logarithm of arbitrary units in projected density. 

ter, while there are Nx, = 171019 and iVhr = 729000 DM 
particles m total, with masses Mir, dm = 2.61 X lO^Uf© 
and Mhr,DM = 3.26 X IQ^^Mq and the same number of 
gas particles with masses Mir, gas = 3.93 x 10^°Mq and 
Mhr.gas = 4.92 X IffMQ. 

3.2 Mapping the SZ effect 

To render the SZ effect we have built a parallel imager ca- 
pable of integrating any physical quantity along the line of 
sight in the simulation. The imager includes the SPH kernel 
of GADGET and tree neighbour searching so that physical 
quantities are well computed at each 3D pixel of the simula- 
tion volume. The tree is filled with the Sph particles and we 
use this to compute the local physical effects with the same 
algorithm that is used to estimate hydrodynamical forces on 
a particle. We then integrate up the results along the line 
of sight (here the Z direction). Figures 3 and 4 show the 
p^tSZ effect and the kpSZ effect respectively. The figures 
show a box 70 Mpc wide and we have projected the full 
thickness of the simulation. Furthermore, to compare with 
the distribution of the projected gas velocities we plot these 
in Figure 5, where we have projected only the gas veloci- 
ties and gas densities from a 3 Mpc thick slice containing 
the highest density region of the simulation i.e. the centre of 
the cluster (the width is again 70 Mpc). We select isoden- 
sity levels in Figure 2 and show them on the other figures 
to make a comparison between the polarisation signal and 
the possible presence of clumps of dense gas. The colour 



scales for the polarisation plots have been kept the same 
to facilitate the comparison. Following Nagai & Kravtsov 
(2003) we have tested the effect of convolving the resulting 
maps with a gaussian kernel corresponding to 1 arcminute 
FWHM. Putting the cluster at z=0.1 this would represent 
a physical scale of 77.5 kpc. In practice, we find that 
convolution leaves the large scale features of the maps un- 
changed but removes some details in the zoomed map of the 
kSZ effect shown in Fig. 6 which has a 20 Mpc width. As a 
result, we present all maps but that of Fig. 6 unconvolved 
for simplicity. 



4 DEPENDENCE ON PHYSICAL 
PARAMETERS 

In this section, we first discuss features of the p^tSZ effect 
associated with resolution issues in the core of the cluster. 
Then we compare the p^tSZ effect to the kpSZ effect. The 
location and intensity of the bright annulus of the p'^tSZ 
effect can provide useful cross-checking information on X-ray 
observations and in particular on the validity of the /3-profile 
as an accurate description of the projected inner density of 
the gas, while the kpSZ effect can constrain the projected 
density-weighted tangential velocity of the gas. 



4.1 Dependence of the p'^tSZ effect on cluster 
internal gas profile 

4.. 1.1 Effect of the internal structure on the intensity of 

the p^tSZ effect 

As we have pointed out above, the computation of the polar- 
isation due to double scattering involves two integrals. The 
first one integrates the effect along the line of sight while the 
other calculates the anisotropy of the incoming radiation at 
each point of the line of sight. In regions of high electron den- 
sity, the two consecutive scatterings will happen at almost 
the same location and we should therefore expect the po- 
larisation term due to double scattering to contain a strong 
dependence on the spatial distribution of electrons in the 
cluster. In fact, if we calculate the Stokes parameters for a 
toy model (here a /3- model for the g£is density profile), we 
find that the overall signal will scale as n^r'^. Prom this ex- 
ample we infer that the double scattering term will be very 
sensitive to the internal gas structure of the cluster like the 
inner radial density slope or the possible presence of a core 
which will affect both and no. If the angular resolution is 
sufficient, the p^tSZ effect could detect clumps of dense cold 
gas orbiting in the ICM such as remnants of minor merg- 
ers or even cold fronts provided their density contrast with 
respect to the surrounding hot ICM is high enough. 

We have already noted in paragraph 2.2.3 that the two 
scattering terms will show a very particular polarisation pat- 
tern with an annulus around the centre of the cluster, and 
that the radius rmax of the max;imum of intensity of this 
annulus is simply related to the core radius if one assumes 
a /3-profile. 
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Figure 3. Map of the p^tSZ effect for the main cluster of the 
highest resolution simulation. The figure shows a box 70 Mpc wide 
and the whole box has been projected. The decimal logarithm of 

the degree of polarisation is colour-coded with the scale provided 
at the top. This scale has been frozen to facilitate the comparison 
with the other maps. Radial lines coming from the central cluster 
are an artifact of the limited number of steps for the integration 
on cp in equation (6) ; they could be avoided but at the expense of 
unreasonable computing time, the CPU scaling as the number of 
iterations. Note the fall-off of the p^tSZ effect as one approaches 
the centre of the cluster: increasing symmetry tends to wash out 
the signal. Note also the residual polarisation on the filament 
leaving to the east: this shows the possibility of the p^tSZ effect 
tracing large-scale structure. 



Figure 4. Maps of the kpSZ effect in the highest resolution sim- 
ulation. The figure shows a box 70 Mpc wide and the whole box 
has been projected. The decimal logarithm of the degree of po- 
larisation is colour-coded with the scale provided at the top (it is 
the same scale as in Fig. 3) . It has been computed using equation 
(4). Contours of the projected gas density have been added to the 
top image to make a precise comparison with previous maps. The 
directions of polarisation due to the kpSZ effect are shown on the 
bottom image. 
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Figure 5. Maps of the projected density and of tlie projected 
gas velocity field in the highest density slice of the simulation. 
The width is 70 Mpc and the thickness 3 Mpc. The length of 
the arrows is proportional to the velocity with a maximum of 
approximately 2000 km s~^ . 



4-. 1.2 Effect of resolution on the estimation of the core 
radius 

As pointed out above, the maximum degree of polarisation 
Pmax of the radiation after the second scattering of the ther- 
mal effect depends significantly on the core radius of the 
/3-profile of the gas. Consequently, Pmax may easily vary by 
a factor 2 for reasonable values of Tc/rvir- This can be advan- 
tageous in the case of multi-frequency observations, where 
one can separate the kpSZ from the p'^tSZ effect and use 
only the latter to constrain the gas density in the central 
region. On the other hand, it shows that it will be difficult 
to disentangle both effects when only a single frequency is 
available. In any case, it is essential to have a robust estimate 
of the order of magnitude of the intensity of the polarisation 
signal which is not plagued by resolution issues. In all the 
above, we have ignored the central softened part of the clus- 
ter as we computed the core radius of the gas according to 
the /3-model because it introduces artificial smoothness in 
the derived cluster gas profile. We have noticed that when 
we skip the inner region as one would do to suppress the 
dependence on the resolution of the simulations, the order 
of magnitude of the p^tSZ effect can significantly increase 
if one uses the "numerical approach" of 2.2.2 to estimate 
it. However, the cluster gas profile parameters r^ and p, 
obtained in this case with the numerical approach are rea- 
sonably close to those resulting from best fitting a /3-model, 
showing that the amplitude we derive for the p^tSZ effect 
inside massive clusters is robust in this case. 

To summarise, even if a very precise computation of 
the intensity of the p^tSZ inside clusters will require accu- 
rate knowledge of the gas profile, our simulations conversely 
prove that, once it is singled out in observations, it can prove 
a powerful tool for constraining the gas profile. 



4.2 Comparison of the two intensities 

We have confirmed numerically the analytical predictions of 
S80 and S99 that finite optical depth effects are not neg- 
ligible and that second scattering of the thermal SZ effect 
may be the dominant source of SZ polarisation other than 
the intrinsic CMB quadrupole in the case where the optical 
depth is high enough. In this paragraph, we compare the in- 
tensities of the p'^tSZ and the kpSZ effects first with respect 
to their spatial distribution in the cluster, then with respect 
to the observing frequency. We also discuss the directions of 
polarisation of the p^tSZ and the kpSZ effects. 



35 
») Mpc 



Figure 6. Map of the kSZ effect in terms of AT/Tcmb- The 
change in temperature is shown with a colour linear between ~ 
— 10~^ and ~ 10~^. Density contours have also been overplotted 
on the map but please note that it is now only 20 Mpc wide. Small 
scale details have been smoothed out as we have the convolved 
the image with a Gaussian beam of FWHM 1' 



4-2.1 Spatial variation 

The peak of the kpSZ effect is located near the centres of gas 
clumps or close to the centres of clusters but seldom directly 
on them. As the kinetic polarisation refiects mostly the part 
of the tangential velocity of the gas which is coherent over 
some scale this is not unexpected because the gas inside 
discrete clumps or in the core region of the cluster may have 
thermal velocities largely in excess of the bulk motions while 
gas accretion on these structures can reach bulk velocities 
of order of 1000 km s~^ which does more than compensate 
for the difference in projected density. 

When the accreting gas reaches the local density max- 
ima its velocity distribution becomes isotropic, hence reduc- 
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ing the signal. This is the case, for instance, for the infaU 
one can observe just under the central cluster (in green), at 
a; ~ 35 Mpc and j/ ~ 35 Mpc in the bottom map of Fig. 4: 
it is associated neither with a local peak of density nor with 
the virial radius of the cluster. It is plain gas accretion, in 
other words, a coherent inflow. On the opposite side of the 
cluster, a clump above the central regions (clearly seen in 
Fig. 2 at a; ~ 45 and y ~ 55 Mpc) is almost absent from the 
kSZ map. If the gas density profile in the very central region 
of a halo is sufficiently symmetric, any similar steep varia- 
tion in the tangential velocity of the gas (accretion or out- 
flow because of AGN or SN feedback) would lead to possible 
detection of the inner ring in the kpSZ maps. Farther out, 
dense cold fronts resulting from recent mergers could also 
make a strong signature in kpSZ maps, as they are thought 
to be contact interfaces separating the hot gas of the ICM 
and the cold gas of the subclumps as the two phases move 
with transsonic velocity with respect to each other. 

It may be noticed that the intensity map for the kpSZ 
effect is sharper than for the p'^tSZ effect. This is mainly 
because the polarisation follows the locations where the gas 
accretes on the high density regions or reaches high veloci- 
ties, and because changes in the projected velocity field of 
the ICM can be very abrupt and combine with variations 
of the projected gas density. On the other hand, the p^tSZ 
effect mostly follows the projected density and temperature 
with the latter being much smoother than the variation of 
velocities. 

4-2.2 Frequency variation 

The different frequency-dependence of the p^tSZ and kpSZ 
effects will introduce changes in the degree and direction of 
polarisation of the resulting combined signal as the observ- 
ing frequency changes. 

As an example, wo briefly discuss the effects expected 
as the observation frequency changes from Xi, < 3.83 to 
Xv > 3.83, looking at infalling gas on the outskirts (r <^ rvir) 
of a massive cluster. 

For Xi, < 3.83 the direction of polarisation of the p^tSZ 
effect is approximately orthogonal to the direction of the 
kpSZ effect. As the cluster is the geometrical centre for both 
the accretion and the tSZ effect, it can be viewed as a point 
source far from the core. Therefore both effects share the 
same symmetry but one is rotated 90° with respect to the 
other. Because they have similar intensities they may coun- 
terbalance one another and the resulting degree of polari- 
sation can vanish. Conversely, if > 3.83 the two effects 
are coherent and oriented in the same way so they will am- 
plify: the signal can be twice as large. (Recall that in all our 
figures we have assumed a,n Xv = 3.) 

4-2.3 On the direction of the resulting polarisation 

The direction of polarisation of the kpSZ effect is clearly or- 
thogonal to the direction of the tangential velocity of the gas. 
In our simulations the cluster is massive enough to drag the 
gas from even the distant edges of the box: on the velocity 
fields, arrows in the outskirts of the box are pointing towards 
the central cluster. It is interesting to note that along major 
structures like filaments this favoured orientation is also sig- 
nificantly enhanced: coherent motion of the gas is amplified 



over such regions. Conversely, at this level of resolution, the 
orientation of polarisation may point toward a completely 
arbitrary direction over the close outskirts of the cluster: 
the signal generated by gas accretion can interfere with that 
due to the presence of slowly moving, high density clumps 
around the cluster. 



5 OBSERVATIONAL PROSPECTS 

In this section, we go beyond using SZ polarisation to con- 
strain only the dynamics and density profile of the ICM. We 
first discuss how the detection of the kpSZ polarisation could 
be used as a tracer of collapsing objects. Then, we present 
the p^tSZ effect as a means of probing external structures. 
Finally, we discuss the contaminating polarisation induced 
by other dominant SZ effects, by the CMB polarisation, by 
dust, and by the cluster magnetic field. 

5.1 The kinetic polarisation term as a probe of 
collapsing structures 

The key feature of the polarisation due to the kinetic term is 
its dependence on the tangential component of the velocity. 
In the situation where a structure is collapsing under the 
effect of its own gravity, the velocity vectors point toward 
the centre of the structure. In this situation it is possible to 
calculate the integrated effect along the line of sight at the 
angular direction 6 from the centre. The degree of polarisar 
tion scales as: 

p oc nf,{v cos{a))^ dl (13) 

where v cos(a) gives the tangential component of the veloc- 
ity vector. We can substitute cos^(a) = d^ /(d^ + l^) where 
d = 6Ds,ng with Dang the angular diameter distance to the 
cluster. This integral over the line of sight of the tangential 
component times electron density will reach a maximum at 
some point between the centre of the structure and its rar 
dius but this maximum will never happen at the centre of 
the structure since in this case the tangential component of 
the velocity is very small (cos(a) — > 0). 

Although the above discussion was focused on a clus- 
ter, the same discussion will be valid for a structure which 
is still collapsing, such as proto-clusters. This opens the ex- 
citing possibility of detecting the clusters even before their 
collapse time. This is because the collapse polarisation term 
is proportional to the collapse velocity squared times the 
optical depth which is the same for a cluster already formed 
and for its progenitor in the collapsing process, and also be- 
cause temperature does not play any role. In fact, we only 
require that the gas cloud is ionised and collapsing. It could 
be possible to see this gas in SZ polarisation forming a ring 
around the centre of the cloud before it starts to emit X- 
rays or before its temperature is high enough to create a 
substantial thermal SZ effect. At distances r > r^, the col- 
lapsing term will dominate the effect due to the bulk motion 
of the cluster. This comes from the higher collapsing veloc- 
ities expected for a typical cluster when compared with its 
expected bulk motion. An example is shown in Fig. 7, where 
the shape of the kpSZ effect expected for a spherical proto- 
cluster collapsing with a homogeneous collapsing velocity of 
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Figure 7. The solid line shows the polcirisation due to the col- 
lapse of a cluster (rc = 150 kpc and Ue = 10~^gr/cm^) with 
a homogeneous collapsing velocity of 800 km/s. The dotted line 
shows the effect due to the bulk motion of the cluster with re- 
spect to the Hubble flow for a tangential velocity of 300 km/s. 
The real polarisation amplitude is the product of the plot times 
the corresponding frequency dependence. 



800 km/s (solid line) is compared to the effect of bulk mo- 
tion of the cluster assuming a tangential peculiar velocity of 
300 km/s (dotted hne). 



5.2 The finite optical depth effect as a 
polarisation generator 

S80 have noted that CMB polarisation due to double scat- 
tering effects can bring out gas structures much less over- 
dense than the virialised regions of massive clusters: smaller 
groups and filaments with sufficient optical depth may be- 
come sources of polarised sub-mm radiation, although to a 
very small degree. This is because the second scattering of 
the primary CMB modified by either the tSZ or the kSZ 
effect does not depend on the local temperature or on the 
local velocity, but only on the local optical depth. As a re- 
sult, if it is possible to measure the few nK of the p^tSZ 
effect up to the outskirts (e.g., rvir) of massive clusters, it 
will be feasible to detect filaments as clusters would provide 
the tSZ input to generate the additional CMB anisotropy. 
An example of this effect is the left filament leaving the 
cluster in the NW direction in Fig. 3 (see at a; ~ 20 Mpc 
and ?/ ~ 40 Mpc), which is less apparent in the kpSZ map 
(Fig. 4). This is in contrast to the mapping potential of the 
tSZ effect in temperature which is essentially restricted to 
lines of sight where the integrated Compton parameter y 
reaches substantial values because at each point, the high 
local temperature must contribute to significant pressure. 

To summarize, it is possible to detect filaments, first, 
by mapping an area of the sky and second, by modelling the 
distribution of nearby clusters: it may be possible to pre- 
dict the direction of polarisation expected on a given region 
of the filament. Here, simulations including photoionisation 
from the UV background that strongly affects the IGM are 
necessary to make a quantitative assessment. 



5.3 Contamination effects 

The kpSZ and p^tSZ effects are probes that follow respec- 
tively the dynamics and the inner density profile of the ICM 
but it is not clear whether they may be detectable in the 
near future. Intrinsic CMB polarisation, SZ effects (scatter- 
ing of the intrinsic quadrupole, and possibly even octupole: 
Challinor et al. 2000), the tangential kinetic effect, and sec- 
ond scattering of the tSZ and kSZ effects), foregrounds such 
as Galactic dust and Galactic synchrotron emission or radio 
point sources with synchrotron emission all contribute to the 
polarisation of the incoming sub-mm signal. These effects 
have a different frequency dependence except for scattering 
of the CMB intrinsic quadrupole which depends linearly on 
in the same way as does second scattering of the kSZ ef- 
fect. The obvious way to separate components will be with 
multi-frequency observations. 

However, the situation could be tractable even at a 
single frequency, provided one obtains an image with suffi- 
ciently high resolution to avoid point sources. Galactic dust 
and Galactic synchrotron emission induce a high degree of 
polarisation, but dust maps will enable one to correct for 
some of the contamination or at least to set a lower galactic 
latitude for the observations. Magnetic fields in the Galactic 
ISM will Faraday rotate the plane of polarisation, however, 
Faraday rotation depends strongly on the frequency. Be- 
cause our present knowledge of the Galactic magnetic field 
is still crude this constraint may also translate into a mini- 
rimm Galactic latitude for observations or to a lower bound 
for the operating frequency. 

Removal of the intrinsic CMB polarisation is feasible, as 
their fluctuations are typically on scales larger than the ones 
considered here (i.e. greater than a few arcminutes). An ap- 
propriate filtering of the maps can safely remove most of this 
signal while keeping the small scale polarisation in clusters. 
At low redshift when the quadrupole remains close to the 
one observed at z=0, the polarisation component in galaxy 
clusters due to the intrinsic CMB quadrupole will vanish in 
four directions in the sky which have been measured with 
some accuracy by WMAP. Furthermore, a measurement of 
the CMB quadrupole would in principle enable a modelling 
of its induced polarisation signal in galaxy clusters which 
could be removed from the data. 

Note that although complex large scale or intracluster 
magnetic fields may again Faraday rotate the plane of po- 
larisation, the rotation angle at frequencies of ~ 300 GHz 
is negligible: according to Ohno et al. (2003), the angu- 
lar rotation of the CMB quadrupole is typically of order 
Arm ~ 1 - 10°{WGHz/i'f on the hne of sight through a 
massive cluster. 

Another possible source for generating the CMB 
anisotropy necessary to obtain polarisation inside clusters 
is through gravitational lensing effects (Gibilisco 1997). For 
the fairly massive Abell cluster A576 (Mvir ~ 1.110^^ 
h''^ Mq according to Girardi et al. 1998), (Gibilisco 1997) 
finds that the expected contribution from gravitational lens- 
ing to the CMB polarization is more than an order of mag- 
nitude less than the expected contribution from the p^tSZ 
effect. Note that (Gibilisco 1997) also considered two other 
gravitational effects due to the gravitational bound or to the 
possible gravitational contraction or expansion of the clus- 
ter, and found that for A576 both effects can be two orders 
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of magnitude greater than the p'^tSZ effect. Ifowever, these 
estimations have been obtained in the frame of a "two-steps" 
vacuole model (Nottale 1984) which is not supported by sim- 
ulations of large-scale structure. This goes beyond the scope 
of the present paper, but it scorns necessary to quantify the 
impact of the last two gravitational effects more precisely in 
future work using simulations. 

With multi-frequency observations and extensive con- 
tamination removal, we conclude that the kpSZ effect can be 
used to inoasurc the projected, density-weighted tangential 
velocity of the ICM in massive clusters, to the level shown 
in Fig. 4. On the other hand, the kSZ effect can be used 
to measure the line of sight component of the projected, 
density-weighted radial velocity of the gas. An example of 
the kSZ effect is shown in Fig. 6, which is a 20 Mpc wide, 10 
kpc thick slice of the central cluster. It is obviously possible 
to compute the velocity of major clumps on the map. The 
issue is whether polarimeters can reach a sensitivity of order 
of a few nK. 



6 CONCLUSIONS 

We have carried out high resolution hydrodynamical simu- 
lations of the formation of a massive, Perseus-size cluster to 
check whether the polarisation of the CMB induced by the 
SZ effect at z ~ can provide information about the dynam- 
ics and inner density profile of the ICM. The contribution 
to the observed polarisation due to the SZ effect can be 
split into 4 major contributions: scattering of the intrinsic 
quadrupole, scattering of the velocity-induced quadrupole 
and second scattering of respectively, the thermal and ki- 
netic SZ effect. The different frequency dependences will 
enable them to be separated using multi-frequency obser- 
vations. We have not discussed here the polarisation due to 
the CMB quadrupole as it is straightforward to compute, 
only depending on the projected electron density, and can 
also be avoided in four directions of the sky, although we 
note that it can provide cross-checking information on the 
gas distribution inside the cluster. We have focused on (1) 
the polarisation due to the velocity-induced quadrupole and 
on (2) the polarisation due to second scattering of the tSZ 
effect because they are possible strong constraints on respec- 
tively the projected density-weighted tangential velocity of 
the gas and on the projected density of the gas within the 
cluster core radius, and we have shown that our numerical 
simulations can make accurate estimates of both amplitudes. 

When the kpSZ effect is combined with the kSZ effect, 
which probes the projected, density-weighted radial veloc- 
ity of the gas, it is possible to deduce information about 
the 3D dynamics of the ICM, at least in regions where the 
projection along the line of sight washes out only a limited 
fraction of the signal. This is the case, for instance, at the 
projected location of dense knots of cold gas orbiting the 
ICM and which reside inside the virial radius as a result of 
past mergers. We plan to tackle how well a typical 3D veloc- 
ity field can be recovered from both the projected kSZ and 
kpSZ information in a subsequent paper. 

Interestingly, the kpSZ effect can also prove very useful 
for mapping collapsing structures, where the tangential com- 
ponent of the collapse velocity can reach values higher than 
the bulk velocity induced by large-scale structures. This is 



further facilitated as the kpSZ term does not depend on the 
temperature of the gas and the collapsing gas is typically 
not hot enough to be studied with X-ray observatories un- 
til structures are well advanced in their formation process. 
Note also that the optical depth of a collapsing structure is 
redshift independent. As a result, the kpSZ effect allows to 
map these proto-clusters higher in redshift than do other, 
temperature dependent effects. 

Although the computation of the polarisation (2) due to 
second scattering of the tSZ effect is intricate, wc have con- 
firmed with more detailed modelling, and for the first time 
in a realistic cosmological setting, that this effect can reach 
the level expected from the scattering of the quadrupole in- 
duced by the tangential motion of the gas. This p'^tSZ effect 
scatters the first anisotropies created by the thermal SZ ef- 
fect occurring inside massive clusters, and because the tSZ 
anisotropy can be seen much further out without diminu- 
tion, the resulting polarisation signal can trace relatively 
well structures like groups, filaments and the cosmic web 
which are much cooler and also less ovcrdcnsc. 

The challenge is now for CMB polarisation-sensitive in- 
struments to reach the 10 nK level typically predicted for 
the SZ effect towards massive clusters. 

High resolution versions of the maps shown are available 

from the website: 

http://www-astro.physics.ox.ac.uk/~gxl/sz_images/ 
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